Quantum localization and cantori in chaotic billiards 
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We study the quantum behaviour of the stadium billiard. 
We discuss how the interplay between quantum localization 
and the rich structure of the classical phase space influences 
the quantum dynamics. The analysis of this model leads to 
new insight in the understanding of quantum properties of 
classically chaotic systems. 

PACS number: 05.45. +b 

The study of quantum mechanics of complex systems 
at the light of the chaotic behaviour of the corresponding 
classical systems has greatly improved our understanding 
of quantum motion M. For example, the possibility put 
forward long ago Q that Random Matrix Theory may be 
a convenient tool to describe spectral properties of clas- 
sically chaotic systems rests now on more solid grounds 
HQ. However, in spite of the progress in recent years and 
the growing interest in the so-called "quantum chaos" we 
are still far from a satisfactory understanding, as a great 
variety and, to some extent, unexpected rich behaviour 
of quantum motion continues to emerge j^] for which a 
satisfactory explanation is required. For example, the 
phenomenon of quantum dynamical localization discov- 
ered 20 years ago in systems under external perturbations 
Q and now experimentally confirmed |Q, mainly rests 
on numerical computations and on qualitative consider- 
ations. Only few mathematical results exist and it is not 
clear whether existing semiclassical theories can account 
for this important feature. Even less understood is the 
mechanism of dynamical localization in conservative sys- 
tems Js|. Billiards are very convenient models to study 
since they can display a rich variety of dynamical behav- 
ior from completely integrable (e.g. the circle) to weakly 
mixing (e.g. the triangle with irrational angles) j^], to 
completely chaotic with power law decay of correlations 
(the stadium), up to exponential decay of correlations 
(dispersive billiards) ]l0| . In addition they can studied 
numerically with great efficiency and high accuracy (here 
we are able to compute accurate eigenvalues and eigen- 
functions of the desymmetrized billiard with sequential 
number up to 10 7 ). They can be, to some extent, stud- 
ied analytically and in laboratory experiments |pd[-|l3[| . 
Also they may be relevant for technological applications 
such as the design of novel microlasers or other optical 
devices fpjj . 



Recently localization has been shown to take place in 
the stadium billiard p"5| and on other similar models 
||i~G| fiol . It is associated to the fact that for small pertur- 
bations of the circle, the angular momentum undergoes a 
classical diffusive process and quantum effects may lead 
to suppression of this diffusive excitation. 

In this paper we study the localization phenomenon 
and the structure of eigenfunctions (EF) as one moves 
from the perturbative regime to the ergodic, delocalized 
regime. The rich variety of classical phase space deter- 
mines a quite complicated quantum structure. Indeed the 
classical motion in the stadium billiard can be described 
by a discontinuous map of the saw-tooth type. This map 
is known to have cantori |pof which may act as barriers 
to quantum motion plfl . This effect has been discussed 
in [p2| and recently confirmed in numerical computations 
on the saw-tooth map on the cylinder jl8| . 

In the following we discuss how the combined presence 
of the cantori structure and of quantum dynamical local- 
ization acts on EFs until the regime of quantum ergodic 
behaviour is reached. We consider the motion of a free 
point particle of unit mass and velocity v (energy E = 
v 2 /2) bouncing elastically inside a stadium-shaped well: 
two semicircles of radius 1 connected by two straight line 
segments of length 2e. The classical motion, for arbitrary 
small e, is ergodic, mixing and exponentially unstable 
with Lyapunov exponent A ~ e 1 / 2 . It can be approxi- 
mated (up to 0(e)) by the discontinuous stadium-map 
]l5| L n+1 = L n - 2esin6>„sgn(cos6'„) v /I - Ll,6 n+1 = 
On + — 2 arcsinL n+1 for the rescaled angular momen- 
tum L = lj\/2E where I = r A v, and the polar angle 
(identical to arc- length for small e). From rigorous 
results on the sawtooth-map Q and from the stadium 
dynamics it can be shown that the angular momentum 
(for small e) undergoes a normal diffusive process with 
diffusion rate 

D = ((l n - lo) 2 )/n\ n>>1 « 2e b ' 2 (2E - I 2 ) (1) 

(Notice that diffusion rate D depends on the local value 
of angular momentum). The power 5/2 in eq. (0) is due 
to the existence of cantori which form strong obstacles 
to phase space transport. In ref [l5| ] the phenomenon of 
quantum localization has been shown to take place in the 
stadium billiard leading to strong deviations from RMT 
predictions. However the dependence of the localization 
length on system parameters is not known and in partic- 



ular we do not know if, and to what extent, the presence 
of cantori will influence the quantum dynamics. Indeed 
it has been conjectured |^2| that cantori act as perfect 
barriers for quantum motion provided the flux through 
cantori is smaller than a Planck's cell, T < 2ttT%. On the 
basis of results on the sawtooth-map [^(j) we can estimate 
the flux T — phase space area transported through can- 
tori per iteration (bounce with the boundary) — which 
is here independent of the winding number of the reso- 
nance, and, for small e, it is given by T ps (2_E) 1 / 2 e 3 / 2 
which leads to the cantori border 

e c = k~ 2 / 3 , (2) 

where k = \J2E is the wavenumber. For e < e c , 
[x := fee 3 / 2 < 1), cantori act as perfect barriers, and 
the quantum system looks as if classically integrable. It is 
therefore expected that the localization length t of eigen- 
states in angular momentum variable I must be of the 
order of the size of cantori. This size, in rescaled angular 
momentum variables, averaged over all the resonances, 
can be estimated from the exact results on sawtooth map 
p3j, namely p — ce, where c ps 12 is a numerical con- 
stant, (c = 10 for e = 0.05, and c = 15 for e = 0.005). 
The fact that c slowly increases with decreasing e is due 
to the presence of the cantorus along the separatrix of 
2:1 resonance (around L = 0) which has a larger size, 
p(2, 1) ps y/i. In Fig|l] we show the classical structure 
of cantori (e = 0.003) in phase space around the largest 
2:1 resonance and the associated quantum eigenstate. In 
this regime the (average) rescaled localization length of 
eigenstates a = £/£ max — ijk is indeed found to be equal 
to the (average) size of cantori (see Fig.|],FigJ|), er = p. 
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FIG. 1. Time evolution of a single classical orbit, followed 
up to 20, 000 bounces, for the classical billiard with e = 0.003. 
The orbit is initially started in the middle of the largest 'is- 
land' ( L = 0,0 = 7r/4 - 0.0016) (a). Angular momentum 
probability distribution Pk(l) of the corresponding eigenstate 
* fc with e = 0.003 and eigenvalue k = 5999.8166. As it is 
seen the state is uniformly distributed over the cantorus in 
the main island (b). 

For e > e c , (x — fee 3 / 2 > 1), when the flux trough turn- 



stiles becomes larger than one quantum 2irh, the cantori 
do not act any more as barriers for quantum dynamics 
and the quantum motion starts to follow the classical 
diffusive behaviour up to the quantum relaxation time 
tn (break time), which is proportional to the density of 
operative eigenstates. For t > t R instead, the quantum 
dynamics enters an oscillatory regime around the station- 
ary localized state with a localization length I. Therefore, 
t R = adN/dE : where a = £/£ max = i/k, and N = E/8 
(from Weyl formula). From the diffusion law (1) we have 
I 2 ps Dt R /T = o-D/T, where T ps E^ 1 / 2 is the average 
time between bounces. Thus we obtain a simple expres- 
sion for the rescaled averaged localization length a 

a ps D/k ps ae 5/2 k = aex. (3) 

where a ps 1.7 is a numerical constant. However, to the 
above expression, we need to add the average size p of 
cantori. Therefore, for £ > 1, the actual expression for 
the localization length will be given by 

cr = p+ (1 -p)ae(x - 1) (4) 

which takes into account also the fact that we need to 
rescale the total size of angular momentum space, and 
that for x = 1, a = p. Eigenstates become delocalized 
(ergodic) when a = 1, giving the ergodicity border e = e e , 

e e ps {ak)- 2 ' b (5) 

which agrees with the results of ]l5| . The cantori border 
can actually be observed if it is below the ergodic border 
and above the perturbative border e p . The perturbative 
border e p is given by the condition that e should be large 
enough, e > e p , to couple two neighboring eigenvalues of 
angular momentum, which is equivalent to the intuitive 
condition of comparing the deformation e and the De 
Broglie wavelength, which leads to 

e P = k~\ (6) 

Therefore, for sufficiently large k, we have, e p < e c < e e . 
In this situation it is natural to expect that cantori will 
influence the localization process and we may have here a 
nice possibility to study the effect of cantori in quantum 
mechanics. 

In order to check the above predictions we numer- 
ically computed quantum eigenfunctions ^(r) of the 
stadium billiard (solutions of the Schrodinger equation 
(V 2 + k 2 )^k = where % = 1) by expanding them 
in terms of circular waves (here we consider only odd- 
odd states) ^>k(r) = J2tLi a sJ2s(kr)sm(2s9). Eigenval- 
ues k = k n and the associated coefficients a s have been 
computed very efficiently p9| by minimizing a special 
quadratic form defined along the boundary of the bil- 
liard J24| . The coefficient a s is proportional to the prob- 
ability amplitude of finding angular momentum equal to 
I = 2s ,PkQ) = \{l - 2s\V k )\ 2 = \a s \ 2 drr\J 2s {kr)\ 2 cx 

KlVi-p/fc 2 . 



Quantum localization in the stadium is not exponen- 
tial like for kicked rotor or smooth diffusive billiards jlj) , 
but rather of power-law type. More precisely, the tails 
of eigenfuctions have been found [[lTJ to decay on av- 
erage, as p(l) ~ \l — |~ 4 . Therefore, in order to 
characterize the localization length of quantum eigen- 
functions we choose the 99% probability localization 
length a rather than the more common inverse partici- 
pation ratio, information entropy or uncertainty of an- 
gular momentum; indeed, the former is quite indepen- 
dent of the nature of tails of the distribution p(l) unlike 
the latter quantities. More precisely, we define a% as 
the minimal number of angular momentum states which 
are needed to support the 99% probability of an eigen- 
state * fe , a k = mm{#A,J2i EA Pk( l ) > 0.99}/(0.725fc). 
The normalization factor has been chosen in such a way 
that crfc = 1 for completely delocalized (GOE) states. In 
Fig. 2 and Fig. 3 we show the dependence of the averaged 
rescaled localization length a = (ak) (averaged over a set 
of consecutive eigenstates ^k with i on a narrow inter- 
val) on the parameter e and on the wavenumber k (up 
to k = 12000, N w 10 7 ). Numerical data agree with 
theoretical predictions (^) with a — 1.7. 

The exact stadium-eigenstates *&k„ may be ex- 
panded in terms of unperturbed quarter-cirlce eigenstates 
$ sm (r) = J2s{k° sm r) sin(2s6 l ) where kP sm are the eigenval- 
ues of the integrable quarter-circle - the zeros of the even- 



order Bessel functions. The matrix i 



= (#. 



kn. 



can 



be easily computed from the coefficients ai |16|. 

It is important to note that unlike for Wigner band 
random matrices || the matrix c" m , (ordered, as usual, 
with increasing quantum wavenumbers k n ,k® m ), has a 
symmetric appearance |]l9| . It has a band structure with 
the bandwidth b ~ k = 5E, independent of e for e < e p , 
where SE is the width of the energy shell. However, be- 
low the ergodicity border, the matrix c™„ is uniformly 
sparse, both in horizontal and vertical directions. The 
effective number of nonzero elements in each row (or 
column) is (on average) equal to the localization length 
£ = ak. In order to illustrate the general structure we 
show in Fig. 4 the probability distributions of a typical 
circle state in terms of eigenstates of the stadium (one 



column of the matrix 



■) and a typical eigenstate 



of the stadium in terms of circle states (one row of the 
matrix). These distributions are strongly sparse and ex- 
tended over the whole energy shell of SE ~ b = k states. 
In Fig. 4c we show, for the same stadium eigenstate, the 
probability distribution p(l) in angular momentum: as 
expected it is strongly localized and non-sparse. The 
above structure is found in the regime of dynamical lo- 
calization and in the regime of cantori localization where 
the quantum system behaves as if classically integrable 
pBf . Furthermore, a similar structure is found also for the 
nearly-circular but chaotic rough billiards with smooth 
boundaries introduced in Jl|] (above the 'Breit-Wi gner 
regime'), where the situation is even more complicated 
due to KAM structure of classical phase space. 

As the parameter e is increased up to the ergodicity 



border e e , sparsity decreases and the quantum angular 
momentum distribution p(l) approaches (apart from fluc- 
tuations) the classical steady-state microcanonical distri- 
bution, p e (l) cx \J2E — I 2 . Notice that the scaling param- 
eter a also controls the deviations from RMT predictions 
of the statistical properties of eigenvalues and cigenfunc- 
tions. 
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FIG. 2. Rescaled localization length a versus the scaling 
variable x = e 3/2 k for five values of e (60 < k < 12,000). 
Each point is obtained by averaging over a large number v of 
consecutive eigenstates (y — 100 for small k and v — 1,000 
for large k). The numerical data clearly show the cantori 
border x = 1. In the cantori region a is constant as expected, 
while for x > 1 the numerical data agree with the theoretical 
prediction (Q) (dotted curves). For large x, the value of a 
approaches the maximal ergodic value a = 1. 
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FIG. 3. Rescaled localization length a versus e for two 
different values of k. The full lines give the classical estimate 
for the average width of cantori. It is seen that below the 
cantori border x = 1 (e = 0.003 for k = 6, 000 and e = 0.0063 
for k = 2, 000) a is proportional to e and independent on k. 
Above the border x = 1 instead, the numerical data follow 
the theoretical estimate (U) (dotted curves). 

In this paper we have discussed a dynamical model, the 
stadium billiard, for which the classical motion is com- 
pletely chaotic (without any island of stability) for any 
value of the control parameter e and we have shown that 
the quantum dynamics instead, exhibits a rich structure 
and different regimes of motion as a function of e and en- 
ergy E. It has been shown that the presence of cantori in 



classical phase space may have strong effect on the quan- 
tum dynamics and leads to a new border which is differ- 
ent from the perturbative and the ergodic border. In the 
regime of quantum cantori (where the phase space flux 
through cantori is less than one quantum) the rescaled 
localization length a = t/k does not depend on energy or 
wavenumber k = y/2E. However, above the cantori bor- 
der, quantum dynamical localization takes place and the 
localization length I is found to be proportional to the 
rate D of classical diffusion in angular momentum. The 
mechanism of localization is strongly connected to the 
sparsity of EFs when expanded on the basis of (unper- 
turbed) circle states (and vice versa). We suggest that 
the above features are typical of the quantum dynamics 
of classically chaotic conservative systems. 

Discussions with R.Prange, R.S.MacKay, J.Keating, 
G. Tanner and I.Dana are gratefully acknowledged. T.P. 
acknowledges financial support from the Ministry of Sci- 
ence and Technology of R Slovenia. 
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FIG. 4. The structure of high-lying eigenstates of the sta- 
dium fc n w 2000 for e — 0.01 in terms of circle states and 
vice versa. In (a) we show the probability distribution of a 
typical circle state (I = 2s = 848, k° sm = 1999.99349) versus 
the wavenumber (energy) of the stadium eigenstates (a col- 
umn of the matrix |c™ m | 2 with fixed s, m), in (b) we show the 
probability distribution of a typical eigenstate of the stadium 
(k n — 1999.91397) versus the wavenumber of circle states 
(a row of the matrix |c" m | 2 with fixed n). In (c) we show 
the same state as in (b) but in angular momentum quantum 
number I (the sum over radial quantum number ^ | c" m | 2 ) . 
Notice that the rescaled localization length here is Ok n = 0.10 
and, correspondingly, only a fraction at n of the total number 
of ~ k states inside the energy shell is actually significantly 
excited (a,b). The dotted curve in (c) gives the ergodic dis- 
tribution p e (i) = (8Tv/k 2 )Vk 2 - P. 
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